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Abstract. We report on a detailed study of the X-ray spectrum of the nearby eclipsing spectroscopic binary YYGem. 
Observations were obtained simultaneously with both large X-ray observatories, XMM-Newton and Chandra. We compare 
the high-resolution spectra acquired with the Reflection Grating Spectrometer onboard XMM-Newton and with the Low Energy 
Transmission Grating Spectrometer onboard Chandra, and evidence in direct comparison the good performance of both in- 
struments in terms of wavelength and flux calibration. The strongest lines in the X-ray spectrum of YY Gem are from oxygen. 
Oxygen line ratios indicate the presence of a low-temperature component (1 — 4MK) with density n c < 2 10 10 cm -3 . The 
X-ray lightcurve reveals two flares and a dip corresponding to the secondary eclipse. An increase of the density during phases 
of high activity is suggested from time-resolved spectroscopy. Time-resolved global fitting of the European Photon Imaging 
Camera CCD spectrum traces the evolution of temperature and emission measure during the flares. These medium-resolution 
spectra show that temperatures > 10 7 K are relevant in the corona of YYGem although not as dominant as the lower temper- 
atures represented by the strongest lines in the high-resolution spectrum. Magnetic loops with length on the order of 10 9 cm, 
i.e., about 5 % of the radius of each star, are inferred from a comparison with a one-dimensional hydrodynamic model. This 
suggests that the flares did not erupt in the (presumably more extended) inter-binary magnetosphere but are related to one of 
the components of the binary. 
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1. Introduction 

The Castor system comprises three visual stars, all of which are 
spectroscopic binaries. The optically faintest of the three com- 
ponents, YY Gem (= Castor C), is of major importance to stel- 
lar evolut ion studies: bei ng an eclipsing spectroscopic binary 
(i ~ 86°; Pettersen 1976[ ) very near or on the main sequence it 
allows for the determination of masses and radii of both com- 
ponents, and for a test of evolutionary tracks. Both stars in the 
YYGem system are of nearly equal spectral type, dMle, in a 
synchronous orbit with a period of 0.8 1 d. YY Gem was the first 
late-type star on which periodic photometric variability was de- 
tected (Kron 1952). frjopp & Evans (1973) have explained sim- 
ilar brightness variations on the flare star BY Dra by rotation of 
star spots, and stars displaying this phenomenon are since then 
termed BY Dra variables. Doppler images of YY Gem have re- 
vealed spots at mid-latitudes on both stars (Hatzes 1995). 

X-rays from the Castor sextuplet were first recorded by 
the Einstein satellite (IVaiana et al. 19811, taillault 1982L Golub 



et al. 1983). Since then the system was observed by virtually all 
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X-ray observatories; see e.g., Pallavicini et al. (1990), Gotthelf 
et al. (1994),|Schmitt et al. (1994)^|Gudel et al. (2001a)|. As the 



spatial resolution of the instruments improved, X-ray emission 
from individual members of the Castor system could be iden- 
tified. At present, all three visual binaries in the sextuplet are 
known to be X-ray sources. The X-ray emission from Castor A 
and B is commonly attributed to their late-type companions be- 
cause the primaries, being early A-type stars without a convec- 
tive envelope and consequently with no ability to support the 
dynamo mechanism, are not expected to be X-ray sources. 

In nearly all of the previous X-ray studies of YY Gem flare 
activity was reported. The luminosity of YY Gem in the soft X- 
ray band ranges between ~ (2—8) 10 29 erg/s. Flares have been 
observed also in other parts of its electromagnetic spectrum. 



Doyle & Mathioudakis (1990) found that the flare activity on 



YY Gem is more than an order of magnitude larger during out- 
of-eclipse times as compared to times when one of the stars 
is eclipsed. On the basis of these U band observations it was 
suggested that the amount of magnetic energy is largest in the 
inter-binary space, leading to frequent energy release in this 
region. A periodicity of the optical flaring rate was reported 
by poyle et al. (1990^ , but has not been confirmed so far. The 



2 



B. Stelzer et al.: X-ray spectroscopy of YYGem 



chromosphere and transition region was examined by Haisch 
et al. (1990), and moderate variability throughout the orbital 
cycle was detected in the emission lines. 

The extraordinary variability of YY Gem has instigated us 
to carry out a coordinated multi-wavelength campaign com- 
prising radio, optical and X-ray observations. As a unique ex- 
ample of an eclipsing spectroscopic binary with nearly identi- 
cal components YY Gem is ideally suited for a study of stel- 
lar coronal structure. In particular we were aiming at looking 
into the possibility of enhanced flare activity in the inter-binary 
space as suggested by Doyle & Mathioudakis ( 1 990) Magnetic 



structures in between the components of binaries have been 
disc ussed also for RS CVn systems (e.g., U chida & Sakurai 
1983). Additional information about the coronal geometry can 
be obtained from observations of X-ray eclipses. The depth and 
duration of the eclipses provide a measure for the spatial ex- 
tent of the corona and contribution of the two components in 
the binary. From a 3D-deconvolution of the X-ray lightcurve 
of YYGem (using an earlier observation by XMM-Newton) 
Giidel et al. (2001a) found that the coronae of both stars in the 
binary are inhomogeneously structured with brighter areas at 
mid-latitudes. Obscuration of coronal features can also be used 
to localize and constrain the emitting region in active stars (see 
e.g., Gunn et al. 1997, 1999, |Schmitt & Favata 1999|). 



As part of our multi-wavelength project the Castor sys- 
tem was observed simultaneously with both XMM-Newton and 
Chandra. While in the EPIC MOS image presented by Giidel 
et al. (2001a) Castor A and B were spatially resolved, our 
Chandra observation allows for the first time to separate the X- 
ray spectra of Castor A and B from each other. We will examine 
the X-ray characteristics of Castor A and B in a related publica- 
tion (Stelzer et al., in prep.). Here we concentrate on the inter- 
mediate and high-resolution X-ray spectrum of YY Gem. We 
emphasize that this observation presents an excellent opportu- 
nity to cross-check the calibration of the dispersive instruments 
onboard both satellites. We analyse the emission line spec- 
trum observed with the Chandra Low-Energy Transmission 
Grating Spectrometer (LETGS) and the Reflection Grating 
Spectrometer (RGS) onboard XMM-Newton independently of 
each other, and derive observed line fluxes for both instru- 
ments. Emission lines in the spectral region of the LETGS and 
RGS include the helium-like triplets from C V to Si XIII, the 
Lyman series of hydrogen-like ions, and numerous iron L-shell 
transitions. By means of time-resolved analysis of spectral pa- 
rameters we check for spectral variability related to changes 
in the activity level of YYGem. We use both 'global' fits to 
the EPIC X-ray spectrum from 0.3 to lOkeV and modeling of 
individual emission lines in the range 2-175 A (resolved by 
the gratings onboard both Chandra and XMM-Newton) to infer 
the time evolution of temperature and density of the emitting 
plasma. 

The XMM-Newton and Chandra observations are intro- 
duced in Sect. ^. In Sect. |] we present results from high- 
resolution spectroscopy, involving a detailed comparison of the 
contemporaneous LETGS and RGS spectrum. Time-resolved 
spectroscopy is discussed in Sect. We draw conclusions on 
the coronal structure by means of loop modeling (Sect. |5|), and 
summarize the results in Sect. ||. 



2. Observations 

YYGem was observed by bot h Chandra (Obs-ID 28) and 
XMM-Newton ( |obs-id 0112880| 801) on Sep 29/30 2000 for a 
total observing time of 59 ksec and 55 ksec, respectively. Both 
observations were part of the Guaranteed Time program. X-ray 
lightcurves obtained with the individual instruments demon- 
strating the times of observation are shown in Fig. [j]. 

We have computed the orbit al phases for the time o f ob- 
servation using the ephemeris of Torres & Ribas (2002) . Two 
large flare events, the secondary eclipse, a flare-like feature af- 
ter secondary eclipse, and an extended "high state" at the end 
of the Chandra observation are seen. The structure of the X-ray 
lightcurve resembles that of the earlie r XMM-Newton obs erva- 
tion of YYGem made in April 2000 ( |Gudel et al. 200 la[ ). The 
larger of the two flares and the feature after secondary eclipse 
are seen at similar orbital phases in both observations. 



2. 1. XMM-Newton 

The XMM-Newton observation on Sep 29/30, 2000 was car- 
ried out in full-frame mode of both EPIC instruments and with 
the thick filter inserted for both pn and MOS to avoid optical 
loading from Castor. The Optical Monitor was in closed po- 
sition due to the optical brightness of the Castor system. For 
a d escription of the X-ray instruments onboard X MM-Newton 
see Jansen et al. (2001), Jen Herder et al. (2001), Striider et al. 



(2001), Turner et al. (2001), and Aschenbach (2002) 



For the analysis of the EPIC spectrum source photons were 
selected from a circular region centered on the position of 
YYGem. An appropriate outer radius for the extraction of 
source photons is the value where the background subtracted 
radial distribution of integrated counts flattens out, in this case 
~ 26". With this choice we make sure that the data are not 
contaminated by contributions from Castor AB (at a separation 
of 74" from YYGem). A background spectrum was extracted 
from a nearby position with the same area as the source ex- 
traction region. The EPIC data were analysed with the XMM- 
Newton Science Analysis Software (SAS) version from May 
2000. We used the pn detector response made available by the 
hardware team in April 200 1 and the MOS responses released 
in February 2001 (see Sect. 4.2.1 for more details). The RGS 
data were analysed with the SAS version 5.2. CCD PI filter- 
ing allows to separate the orders. We concentrate on the first 
order spectrum. The background spectrum is obtained from a 
region offset from the spectrum of YY Gem and Castor AB in 
the cross-dispersion direction. 

2.2. Chandra 

With Chandra we used the LETGS combining the Low-Energy 
Transmission Grating (LETG) with the High Resolution 
Camera (HRC) (see Brinkman et al. 200C| ). The Chandra obser- 
vation began ~ 2.5 h after the start of the pointing with XMM- 
Newton. After the XMM-Newton observation ended Chandra 
continued to observe for ~ 5 h (see Fig. [I]). 
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Fig. 1. Background subtracted X-ray lightcurve of YY Gem on 2000 Sep 29/30. top - Chandra LETG zeroth order with 200 s 
binsize, bott om - XMM-Newton EPIC pn with 100 s binsize. The orbital phase was computed with the ephemeris by Torres & 
Ribas (2002). The time intervals selected for time-resolved spectroscopy are indicated by horizontal lines above the lightcurves. 



The data obtained from the Chandra science center was re- 
processed following the CIAO 2. 1 science threads[] for apply- 
ing the newest degap-corrections to the HRC-S data as well 
as generating new event level 2 files, in which background is 
reduced significantly by a spatially dependent "light pha filter- 
ing" described in the thread. 

The spectra were then extracted from these event files us- 
ing IDL routines. For the source we used the 'bow tie' extrac- 
tion region as described in the Chandra Proposers' Observatory 
Guide[] (POG). Large areas above and below the source extrac- 
tion region were selected for the background, such that the ra- 
tio between source and background area is 7.5 (see Fig. 1 in 



Burwitz et al. 2001 for an illustration). Towards longer wave- 



lengths A > 48.4 A the background regions become wider in or- 
der to maintain the constant ratio of areas between source and 
background spectral bins. 

3. High-resolution Spectroscopy 

In this section we present and compare results from the 
high-resolution dispersive instruments on Chandra and XMM- 
Newton. We took care that the Castor AB components are ex- 



cluded from the background extraction regions. Lists of all 
identified lines are given in Tables [j] and We derive pho- 
ton fluxes and use line ratios for coronal plasma diagnostics. 
In this context the triplets of He-like ions comprising the reso- 
nance line 7-, the intercombination line i, and the forbidden line 
/, are of particular interest as their intensity ratios R = f ji and 
G = (/ + i)/r provide an estimate for the electron density n e 
and electron temperature T c of the respective lines. The temper- 
ature sensitivity of G is due to the collisional excitation rates 
which have different temperature dependence for the resonance 
line compared to the forbidden and intercombination line. The 
G ratio also indicates whether the plasma is dominated by col- 
lisions or photo-ionization: a strong r line corresponds to colli - 
sionally dominated plasmas (G ~ 1), while r is weak in photo- 
ionized environments with G > 4 (see forquet & Dubau 2000 ). 
We caution that a given line does actually form over an ex- 
tended range of temperatures in a non-isothermal plasma, and 
thus the temperature derived from G may be misleading. The 
number of collisions increases with density. As a consequence 
the upper level of the / transition is depopulated and fills the 
upper level of the i transition. This results in a decrease of R 
with increasing n c . Only in the low- and high-density limits is 
R independent of n c . 



http://asc.harvard.edu/ciao/documents-threads.html 
http://asc.harvard.edu/udocs/docs/docs.html 
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Fig. 2. Time-averaged Chandra LETGS spectrum (left and right side added); total exposure time is 58 ksec, and binsize is 0.04 A. 
Only few emission lines are identified above 100 A (108.4 A, 117.2 A, 128.73 A, and 132.9 A; not displayed here). The region 
within the dotted rectangle (10 — 26 A ) is shown in more detail in Fig. |3| 



3.1. Line identifications 



We extracted and analysed the RGS and the LETGS spectra of 
YYGem as described in Sect. 2.1 



and 2.2. The time-averaged 



LETGS spectrum in the wavelength range from - 100 A is 
displayed in Fig. |2|. We have added left (negative) and right 
(positive) side of the spectrum. Above 100 A only few lines 
are seen, all from highly ionized iron. A close-up of the re- 
gion with the strongest emission lines is shown in Fig. ^| to- 
gether with the time-averaged first order RGS spectrum in the 
same wavelength region. We show the spectrum in units of 
cts/s/bin. Since the RGS and the LETGS observations over- 



lap for about 75 % in time, the measured line fluxes should be 
similar. The count rates are higher for RGS than for LETGS 
throughout most of the displayed spectral region, and demon- 
strate directly the higher sensitivity of the RGS, except at wave- 
lengths < 10 A. The full wavelength range of the RGS extends 
from 6 — 35 A. However, the only lines found in the RGS spec- 
trum of YYGem outside the range displayed in Fig. || are the 
silicon lines between 6 — 7 A, the resonance line of the triplet 
of He-like N VI, and transitions from H-like carbon. 

We identified the most prominent lines in both the LETGS 
and the RGS spectrum. The line centers and fluxes were 




Fig. 3. Comparison of the contemporaneous XMM-Newton RGS 1 and RGS 2, and Chandra LETGS spectra of YY Gem in the range A = 10 — 26 A. Straight horizontal lines 
in the RGS spectra represent gaps due to CCD chain failure or individual chip separation. 
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determined with the CORA^ line fitting application (Ness 
et al. 2001). CORA was designed for low count rate spectra, 
based on Poisson statistics and employs a maximum likelihood 
method to determine the best fit parameters. The background 
is treated as a constant in the vicinity of individual lines. For 
the LETGS spectrum we used Gaussians to approximate the 
line profile. To avoid artificial broadening of lines resulting 
from uncertainties in the wavelength calibration, the positive 
and negative sides of the spectrum were analysed separately. 
Then we checked the agreement, and if adequate, we merged 
the positive and negative sides to increase the signal for the 
analysis of line ratios. For RGS Lorentzians represent a bet- 
ter approximation to the instrumental profile. To improve the 
S/N the first order spectra of RGS 1 and 2 were co-added. Both 
RGS 1 and RGS 2 have a ~ 3— 4 A wide gap in the spectrum as 
a consequence of CCD chain failure (10.4— 13. 8A for RGS 1 
andl9.8-24.2A for RGS 2), and further small (~ 0.1 A) gaps 
due to the finite separation of adjacent chips. 

The line parameters derived from LETGS and RGS are 
summarized in Tables [j] and |[ The results from the positive 
and negative side of the LETGS spectrum are listed separately 
to examine the accuracy of the wavelength calibration. We find 
that the line centers are in excellent agreement. Apart from the 
observed line center and width Tables [j] and || present the to- 
tal number of counts /, the photon flux, and the identification 
and predicted line positi ons adopted from|Mewe et al. (1985' , 
Mewe et al. (1995)| , and |Phillips et al. (1999)| . We used the ef- 
fective areas of the respective instrument to compute the photon 
flux. The photon flux does not depend on instrumental proper- 
ties and can be used for direct comparison of LETGS and RGS. 

In Fig. [| we present a comparison of the LETGS and 
RGS fluxes for individual lines. Our simultaneous observa- 
tions allow to demonstrate that the measurements by these two 
instruments are consistent. Most photon fluxes are in agree- 
ment within one or two sigma. The largest deviations between 
LETGS and RGS occur for the O VII triplet and the Ne X Lya 
line. The deviations in OVII can most likely be attributed to 
an oxygen absorption edge in the RGS which reduces the ef- 
fective area at the position of the O VII triplet, and had not yet 
been corrected for in the res ponse matrix at the time of anal- 
ysis (den Herder et al. 2002). The discrepancy in NeX seems 
to result from a problem in the positive side of the LETGS, as 
its negative side and RGS provide comparable fluxes. Note that 
the photon fluxes in Tables [j] and [2] denote the combined flux 
of both unresolved stars in the YY Gem binary. 



3.2. He-like ions 

The spectral ranges of the LETGS and the RGS comprise the 
He-like triplet transitions of several elements. In the spectrum 
of YY Gem the O VII triplet is by far the strongest and best re- 
solved triplet. The Si XIII triplet is not fully resolved, and the 
Ne IX triplet is blended with iron lines. All other He-like ions 
do not produce sufficient signal. The resonance line of the N VI 




" CORA can be downlo aded from http://www.hs.uni- 
hamburg.de/DE/Ins/Per/Ness/Cora 



RGS flux (10 5 cts/cm 2 /s) 

Fig. 4. Comparison of LETGS and RGS fluxes for individual 
emission lines; 2 a uncertainties. 



triplet is weakly detected in the RGS spectrum, while the inter- 
combination and the forbidden line of N VI are not identified. 
In Fig. U we display the spectral region around the Ne IX and 
the OVII triplets. An iron line on the long wavelength side of 
the NelX forbidden line (at A - 13.83 A) was included in the 
fit. The contamination of the Ne IX triplet stems mainly from 
Fe XIX (at 13.53 A) which is unresolved from the triplet inter- 
combination line. Recall that for both the Ne IX and the O VII 
triplet only one RGS is available. 

The only element with resolved and strong enough lines for 
plasma diagnostics based on an analysis of line ratios is oxy- 
gen. Note that oxygen provides only information about com- 
paratively low temperatures (~ 1 — 4 MK). We determined the 
R and G ratios from the flux in the individual triplet lines. 
Furthermore, we computed the flux ratio between the H-like 
Lya line and the He-like resonance line. This line ratio can be 
used as a temperature diagnostic, because the degree of ioniza- 
tion depends on the temperature. The result is listed in Table |[ 
Fig. H shows that the strength of the r line is comparable to 
the sum of the intensities in the / and the i line. This is typi- 
cal for a collisional plasma. It was shown by Blumenthal et al. 
(1972) that the UV radiation field of the star can modify the 
R ratio in the same way as collisional excitation, namely by 
depopulating the upper level of the / transition to the upper 
levels of the i transition. We made use of the calculations by 
forquet et al. (2001) which include the radiation field of the 
star, to estimate the density n c and formation temperature T of 
the oxygen triplet. However, as YY Gem is composed of two 
M-type stars it is reasonable to neglect the radiation field and 
assume that the relative intensities of / and i reflect the fre- 
quency of collisions and, therefore, represent a direct measure 



Table 1. Results from line fitting of the time averaged first order LETGS and RGS spectra. Positive and negative side of the LETGS are listed separately to demonstrate 
the accuracy of the wavelength calibration. RGS 1 and RGS 2 spectra were combined except for lines which fall in the wavelength ranges corresponding to failed CCDs or 
gaps between chips. Line blends are indicated by quotation marks ("). Lines which are identified but not properly resolved are marked by asterisks (*). For RGS all lines are 
represented by Lorentzian profiles, and for LETGS by Gaussians. The typical error on the line width is ~ 0.005 A. Line identifications are from Mewe et al. (1985), Mewe 
et al. (1995), and Phillips et al. (1999). Lines detected only by LETGS are listed in Table 2. 
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OVIII 


Lya 


18.973 






















28.460 


0.048 


72.6 ± 15.8 


1.6 ± 0.3 


CVI 


Ly/3 


28.446 






















28.776 


0.021 


35.8 ± 10.7 


0.8 ± 0.2 


NVI 


r 


28.787 


33.739 


0.017 


55.0 ± 8.4 


17.5 ± 


2.7 


33.745 


0.026 


83.4 ± 10.3 


23.3 ± 


2.9 


33.732 


0.070 


662.7 ± 31.1 


21.9 ± 1.0 


CVI 


Lya 


33.700 


LETGS +1 order 


LETGS -1 order 


RGS1 




21.614 


0.023 


82.4 ± 10.2 


19.1 ± 


2.4 


21.613 


0.022 


102.8 ± 11.2 


21.8 ± 


2.4 


21.600 


0.047 


475.0 ± 26.0 


14.2 ± 0.8 


OVII 


r 


21.610 


21.812 


0.022 


19.8 ± 6.1 


4.6 ± 


1.4 


21.810 


0.018 


17.6 ± 5.6 


3.8 ± 


1.2 


21.814 


0.045 


99.0 ± 15.3 


3.0 ± 0.5 


OVII 


i 


21.800 


22.115 


0.016 


53.7 ± 8.3 


12.5 ± 


1.9 


22.108 


0.022 


52.7 ± 8.5 


11.5 ± 


1.9 


22.095 


0.044 


332.7 ± 22.3 


10.1 ± 0.7 


OVII 


f 


22.100 


LETGS +1 order 


LETGS -1 order 


RGS 2 




12.150 


0.021 


231.7 ± 16.5 


32.5 ± 


2.3 


12.130 


0.019 


169.8 ± 14.4 


23.1 ± 


2.0 


12.141 


0.044 


736.9 ± 34.0 


21.2 ± 1.0 


Ne X 


Lya 


12.132 


// 


II 






// 


// 


II 


// 




// 


// 


II 


// 


// 


FeXVII 


Ne4C 


12.134 


12.291 


0.019 


33.5 ± 8.1 


4.7 ± 


1.1 


12.266 


0.028 


10.0 ± 7.0 


1.4 ± 


1.0 


12.282 


0.050 


140.1 ± 21.6 


4.0 ± 0.6 


FeXVII 


Ne4D 


12.264 












12.839 


0.027 


22.9 ± 8.1 


3.1 ± 


1.1 


12.836 


0.037 


63.3 ± 15.5 


1.8 ± 0.5 


FeXX 


N16 


12.834 


* 13.467 


0.021 


127.4 ± 13.1 


16.7 ± 


1.7 


*13.465 


0.025 


129.3 ± 13.7 


17.1 ± 


1.8 


*13.451 


0.061 


545.8 ± 35.9 


15.8 ± 1.0 


NelX 


r 


13.448 


* 13.535 


0.024 


40.5 ± 9.4 


5.3 ± 


1.2 


*13.543 


0.030 


59.4 ± 11.1 


7.8 ± 


1.5 


*13.539 


0.072 


248.5 ± 32.8 


7.0 ± 0.9 


FeXIX 


015 


13.524 


// 


// 


// 






// 














// 




NelX 


i 


13.553 


13.720 


0.024 


71.9 ± 10.6 


9.4 ± 


1.4 


13.703 


0.033 


79.1 ± 11.6 


10.3 ± 


1.5 


*13.709 


0.068 


410.1 ± 31.8 


11.3 ± 0.9 


NelX 


f 


13.700 






















*13.831 


0.057 


74.3 ± 21.6 


2.0 ± 0.6 


FeXVII 


Ne3A 


13.826 


17.083 


0.033 


157.4 ± 14.2 


23.4 ± 


2.1 


17.051 


0.035 


152.0 ± 14.1 


20.8 ± 


1.9 


•17.050 


0.053 


419.9 ± 37.5 


12.1 ± 1.1 


FeXVII 


Ne3G 


17.055 


// 


II 










II 






// 


*17.095 


0.052 


393.9 ± 36.8 


11.4 ± 1.1 


FeXVII 


NeM2 


17.100 


24.802 


0.020 


34.3 ± 7.2 


8.0 ± 


1.7 


24.810 


0.027 


57.5 ± 9.0 


12.7 ± 


2.0 


24.788 


0.079 


331.6 ± 23.4 


12.0 ± 0.8 


NVII 


Lya 


24.780 



" Fluxes are not corrected for interstellar absorption. 
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Table 2. Emission lines observed in the time averaged positive and negative LETGS spectra outside the spectral range of or 
not detected by the RGS. All lines are represented by Gaussians. Line identifications are from Mewe et al. (1985), Mewe et al. 
(1995), and|Phillips et al. (1999). 





LETGS positive order 






LETGS negative order 










A 


<T 




Photon Flux 


A 






Photon Flux 




Identification 




[A] 


[A] 


[ctgs] 


[10- 5 SiL^] 


[A] 


[A] 


[cts] 


[10~ 5 -2*-z] 


Ion 


Trans. 


A 




Bern 


s cm 


6.195 


0.027 


ob.z ± 7.o 


3.5 ± 0.7 


6.178 


0.026 


24.7 ± D.o 


2.3 ± 0.6 


Si XIV 


Lya 


6.180 


8.437 


0.019 


18.7 i 6.4 


2.0 ± 0.7 


8.402 


0.019 


20.4 ± 6.D 


2.1 ± 0.7 


Mg XII 


Lya 


8.421 


36.424 


044 


DO o _|_ 7 c 


13.2 ± 2.6 


36.394 


n nos 

U.UZo 


1 ft A 4- K K 


5.6 ± 1.7 


NeX 


Lya(3.) 


36.396 


















S XII 


B6A 


36 398 










36.559 


0.016 


9.8 ± 4.0 


3.0 ± 1.3 


SXII 


B6C 


36.563 


36.912 


0.023 


15.3 ± 4.9 


5.5 ± 1.8 










sx 


N7 


36.900 


37.954 


0.030 


28.0 ± 6.4 


11.2 ± 2.6 


37.934 


0.057 


24.3 ± 6.9 


8.6 ± 2.5 


OVIII 


Lya (2.) 




44.014 


0.020 


30.0 ± 6.3 


4.8 ± 1.0 


44.024 


0.032 


32.5 ± 7.6 


4.4 ± 1.1 


Si XII 


Li6A 


44.020 


44.176 


0.025 


35.1 ± 6.9 


5.6 ± 1.1 


44.169 


0.017 


30.5 ± 6.8 


4.1 ± 0.9 


Si XII 


Li6B 


44.165 


56.920 


0.031 


41.4 ± 7.4 


8.5 ± 1.5 










OVIII 


Lya (3.) 




94.023 


0.059 


70.6 ± 9.9 


28.8 ± 4.0 


93.956 


0.052 


53.9 ± 8.7 


23.3 ± 3.8 


Fe XVIII 


F4A 


93.930 


108.450 


0.058 


42.1 ± 8.2 


22.2 ± 4.3 


108.418 


0.057 


34.0 ± 7.6 


18.6 ± 4.2 


FeXIX 


06a 


108.370 


117.185 


0.058 


46.8 ± 8.6 


26.0 ± 4.8 


117.184 


0.085 


56.6 ± 9.8 


31.9 ± 5.5 


FeXXII 


Bl 1 


117.170 


128.818 


0.049 


33.4 ± 7.3 


33.6 ± 7.4 


128.738 


0.055 


27.7 ± 7.2 


24.8 ± 6.5 


FeXXI 


C6A 


128.730 


132.978 


0.068 


80.8 ± 10.6 


81.7 ± 10.7 


132.952 


0.078 


95.5 ± 11.4 


83.6 ± 10.0 


FeXX 


N6A 


132.850 




13.2 



14.0 




21.0 21.2 21.4 2 



21.8 22.0 22. 
A[A] 



Fig. 5. He-like triplets of Ne IX and O VII measured with LETGS and RGS together with line fits. The continuum close to each 
triplet is approximated by a straight horizontal line. The intensity enhancement to the right of the Ne IX triplet in the RGS 
spectrum is due to FeXVII. The major contamination of the NelX triplet comes from a FeXIX line which is unresolved from 
the triplet intercombination line (see Table [j]). 



of the density. The observed G ratio corresponds to a formation 
temperature between 2 and 3 MK. The formation temperature 
defines the relation between the R ratio and the electron density 
(Porquet et al. 2001). The models place our measurements of 
R (see Table |3|) in the low-density limit, where R is insensitive 
to density, and only an upper limit of n c < 2 10 10 cm -3 can be 
provided. 

The plasma emissivities from the MEKAL code (Mewe 
et al. 1985 and Mewe et al. 1995) applied to the Lya/r ra- 



tio result in a temperature of 3.4 ± 0.1 MK in the LETGS and 
4 . ± . 1 MK in the RGS . These somewhat higher temperatures 



than those derived from the G ratio is naturally explained by the 
higher ionization stage of the involved species. The difference 
between the measurements in these two instruments probably 
arises from the difference in the OVII flux, which is under- 
estimated in the RGS due to calibration problems as outlined 
above. We caution that given these uncertainties the absolute 
values for the line ratios must be treated with caution, but rela- 
tive measurements comparing the same instrument (e.g. as pre- 
sented in Sect. |j) are still meaningful. 
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Table 3. Line ratios R = f /i, G = (/ + i)/r, and Lya/r for 
oxygen, 1 a uncertainties. 





+ 1 


LETGS 

-1 


total 


R 

G 

Lya/r 


2.70 ±0.93 
0.89 ±0.17 
3.73 ± 0.49 


3.01 ± 1.08 
0.70 ±0.13 
3.07 ±0.36 


3.24 ±0.82 
0.77 ±0.10 
3.35 ± 0.29 


RGS 


R 
G 

Lya/r 






3.38 ±0.57 
0.92 ±0.08 
4.98 ± 0.28 



3.3. Optical depth effects 

The general assumption of our plasma diagnostics is that the 
plasma is optically thin, i.e., resonance photons are not scat- 
tered out o f the line of sight a s described by Schrijver et al. 
(1995) and fylewe et al. (1995} To check whether an optically 
thin model is justified we estimate the optical depth by com- 
paring the observed flux ratio of the resonance line of Fe XVII 
at 15.014 A and the adjacent line at 15.265 A of the same ion 
which has a much lower oscillator strength. The same proce- 
dure was carried out for other active stars, e.g. by Mewe et al. 
(2001) and f^ess et al. (2002a) , and they find /15.014//15.265 — 
2.85 ± 0.14 (Capella) and /i 5 .ou//i5.265 = 2.81 ± 0.25 
(Algol). For Proyon we derive a ratio of 2.9 ± 1.7 using the 



RGS measurement given by Raassen et al. (2002). From com- 



parison with laboratory measurements from EBIT (2.8 — 3.2; 



Brown et al. 1998) Mewe et al. (2001) and Ness et al. (2002a) 



concluded that the coronal plasma of Procyon and Algol is op- 
tically thin. From our observation of YY Gem we calculate a 
ratio of 2.8 ± 0.5 for LETGS and 3.2 ± 0.5 for RGS, com- 
patible with the values cited above. This indicates that stellar 
coronae are indeed optically thin, and as the results are simi- 
lar for stars with very different activity levels we conclude that 
optical depth effects generally are negligible in stellar coronae. 
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Fig. 6. Time-evolution of the O VIII H-like Lya line and the 
O VII He-like resonance line, the ratio of these two lines, and 
the R and G ratio of the O VII triplet as observed by XMM- 
Newton RGS. Time intervals as defined in Fig. |IJ 



4. Time-Resolved Spectroscopy 
at High and Low Resolution 

In this section we complement the analysis of the emission 
line spectrum of YY Gem with time-resolved modeling of the 
strongest lines seen with the RGS. The spectra discussed in the 
previous sections comprise the full Chandra and XMM-Newton 
exposure times and represent, therefore, averages over ~ 75 % 
of the orbital period of YY Gem. The spectral shape is, how- 
ever, expected to be variable (i) due to the changing aspect of 
the binary in the course of its orbital motion, and (ii) due to in- 
trinsic variability related to activity on one or both of the stars. 

Next to a time-resolved analysis of individual lines, we 
study the temporal evolution of spectral parameters by a de- 
tailed inspection of the XMM-Newton EPIC spectrum. Giidel 
et al. (2001a) analysed the EPIC spectrum from the earlier 
XMM-Newton observation using MOS 1 and MOS 2 individ- 
ually for 5 different time intervals corresponding to different 
orbital phases of YYGem. Here, we fit the pn and MOS 1 + 



2 spectra obtained in Sep 2000 jointly. This allows to better 
constrain the spectral parameters and perform a more detailed 
systematic analysis of time-resolved X-ray spectroscopy. 

To examine the evolution of the spectrum of YY Gem we 
split the total observing time in several time intervals. In choos- 
ing the time bins we took into account the variability pattern in 
the lightcurve. We found that the LETGS data does not provide 
a strong enough signal for time-resolved analysis, even if left 
and right side of the spectrum are added. Therefore we make 
no further use of LETGS data for time-resolved spectroscopy. 
The selection of the segments for EPIC and RGS can be read 
out of Fig. [j]. 

4. 1. The high-resolution RGS spectrum 

We extracted RGS spectra in five time intervals (see Fig. |l]). 
Segments 2 and 5 are the phases with the highest level of ac- 
tivity. The strongest lines in the RGS wavelength range are ex- 
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amined for all five selected data segments. Some of the weaker 
lines are not identified at all phases. This may be due to vari- 
ability of the source and/or the unequal S/N in the five intervals 
as a result of different exposure times. The only element which 
displays a strong He-like triplet and H-like Lya in all time seg- 
ments is oxygen. For this element we can refine the analysis 
of plasma temperature and density deriving the R and the G 
ratio as a function of time. It must be kept in mind, however, 
that oxygen probes only temperatures between 1 — 4 MK, and 
therefore does not represent the major component of the flare. 

The spectra for the shorter time intervals are characterized 
by poor signal-to-noise. We ensured that the weak intercombi- 
nation line is represented correctly by holding the separation 
between the line centers within the triplet fixed at the nominal 
value during the fitting process. In Fig. |]we show the time evo- 
lution of O VIII and O VII and their line ratios. All ratios are 
characterized by rather large error bars as a result of short expo- 
sure times. A trend towards decreasing R ratio during times of 
enhanced activity (intervals 2 and 5) is evident although only 
marginally significant. To increase the significance of this re- 
sult we merged the three time intervals representing the quies- 
cent state of YYGem, i.e., segments 1, 3, and 4, and obtained 
an R ratio of 3.91 ± 0.79. The difference between this time- 
averaged quiescent R value and the value measured during seg- 
ment 2 is still only a 1 .5 a effect, but is suggestive for a den- 
sity increase as a result either of chromospheric evaporation or 
compression of the coronal material. 

In Fig. [7] we compare the coadded RGS 1 and RGS 2 spec- 
tra during two different time intervals, segment 1 from Fig. [j] 
representing quiescence and segment 5 from the same figure 
representing a phase of strong activity. The bottom panel of the 
figure shows that most of the flare emission is concentrated be- 
tween ~ 8 — 16 A. The low signal does not permit to determine 
whether this emission is true continuum or whether it arises 
from unresolved iron lines which form at high temperatures 
during the flare, e.g., FeXXIV. 

4.2. XMM-Newton EPIC 

Due to its larger effective area the EPIC spectrum can be stud- 
ied at higher time-resolution. We split the EPIC pn and MOS 
lightcurve in a total of 15 time intervals (see Fig. [lj). 

4.2.1 . Photon pile-up 

Photons arriving at the CCD detector can deposit energy in 
more than one pixel. To improve the statistics in the high- 
energy portion of the spectrum we include double events, i.e., 
photons that produced a signal in two adjacent pixels. 

For bright sources (such as YY Gem) more than one photon 
may arrive in one pixel before the detector is read out (photon 
pile-up), and in addition two photons arriving at the same time 
in adjacent pixels are mistakenly identified as a single photon 
distributing its energy on both pixels. Distortions of the spec- 
tral shape due to pile-up and mis-identification of the pixel pat- 
tern distribution can safely be avoided by cutting out the cen- 
tral portion of the point spread function, however, at the ex- 



pense of count rate statistics. In order to compromise between 
correction for these effects and retaining enough events to con- 
strain spectral parameters we examined how model fits change 
if the photons from the innermost part of the source are ex- 
cluded. For different inner radii of the source extraction area 



a 3-temperature VMEKAL model (see Sect. X.2.2) gives sim- 
ilar plasma temperature and abundances within the statistical 
uncertainties. We conclude from this test that the spectral pa- 
rameters are only marginally affected by pile-up in this obser- 
vation, allowing us to use the full source region without having 
to remove the photons from the central area. 

4.2.2. The quiescent spectrum 

Giidel et al. (2001a) described the quiescent X-ray spectrum 
of YYGem by a 3-temperature (3-T) model for thermal emis- 
sion from an optically thin plasma. The existence of a multi- 
temperature corona on YYGem was verified by these au- 
thors also by deriving the emission measure distribution, which 
showed the presence of a broad distribution of temperatures 
ranging from ~ 2 — 15 MK. 

We used a combination of three VMEKAL-models as im- 
plemented in XSPEC 11.0.1 to represent the spectrum in time 
interval t\, when the source was in quiescence. The different 
instruments (pn, MOS 1, and MOS 2) were analysed jointly by 
multiplying the spectral model with a constant normalization 
factor to make up for uncertainties in the absolute calibration 
of the individual detectors. 

We used the EPIC spectrum of phase t\ as baseline for all 
other time segments in which the signal is sometimes lower 
due to shorter exposure time. The model parameters for the 
quiescent spectrum are provided in Table Q The EPIC spectrum 
shows that besides the comparatively low temperatures probed 
by the strongest lines in the high-resolution spectrum, plasma 
in excess of 10 7 K is present. The observed spectrum can only 
be reproduced if the abundances of some elements are free fit 
parameters. The best fit clearly shows subsolar abundances, in 
particular for iron. Low abundances for elements with low First 
Ionization Potential (FIP) were also repor ted from the earlie r 
XMM-Newton observation of YY Gem (see Giidel et al. 2001a). 



4.2.3. Time-evolution of spectral parameters 

At first we checked whether the spectrum of YY Gem can be 
represented by a 3-T model at all times. We assumed that the 
integrated light from the quiescent corona remains unaltered 
during the outburst. Therefore, we held all temperatures and 
abundances fixed and allowed only the intensity to vary, cor- 
responding to changes in the emission measure as could be 
expected due to occultation effects. This way we obtained ac- 
ceptable fits (Xred ~ 1 an d statistically distributed residuals) 
for all time intervals in which the star is 'quiet', i.e. shows 
no strong variability and exhibits intensity similar to the reg- 
ular non-flaring state (£2, tj, t%, tg, iio, and £12). During the 
secondary eclipse (time intervals t$ and tg) the emission mea- 
sures EM2 and EM3 drop because the visible fraction of the 
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Fig. 7. Comparison of the RGS flux spectrum during quiescence (interval 1 from Fig. [T]) and the most active phase (interval 5 
from Fig. |l]). The panel on the bottom shows the difference spectrum with significant excess emission below ~ 16 A. Binsize is 
0.04 A. 



Table 4. Spectral parameters for the quiescent state of YY Gem 
(derived from time interval ti). Normalization constants to 
cross-calibrate the three instruments (pn, MOS 1, MOS 2) are: 
N pn = 1 (fixed), N mosl = 1.0ll°;°l and N mos2 = 1.03±g;°|. 





kT 2 




[kcV] 


n 2i + 05 


64+ 01 


i 7 n+0.30 




EMi 


EM 2 




[io ai cur-'] 


9 24+ 1 - 87 


13.84l°;g 


9 «8+l- 05 




O 


Mg 


Si 




64+0 19 


n 27+°i 2 

-0.07 


fl 47+0-16 

KJAI -0.07 




s 


Fe 


Ni 


Xred (dof) 


n 5n+ 27 

U - OU -0.20 


93+ ' 04 

u - zo -0.03 


u - uu -o.oo 


1.24(658) 



emitting corona has decreased. After the eclipse a strong high- 
energy component (large EM3) is found. These time-intervals 
correspond to an emission feature in the lightcurve (lasting 
from tio to ti2). For t\\, representing the peak of this feature, 
we could not obtain an acceptable 3-T fit. The excess intensity 
in the hottest of the three components is seen already during 
the second half of the eclipse (tg). This could be a hint that the 
emission feature fills in part of the eclipse. 



At all time intervals not discussed in the previous para- 
graph, after fitting the 3-T model to the data a high-energy 
excess stands out in the residuals. This suggests that higher 
temperature plasma is present in addition to the emission from 
the quiescent corona. Adding a 4th VMEKAL component does 
not lead to a significant improvement. Only a 5-T model rep- 
resents an adequate description of the data during the most 
active phases. A 5-T model was also introduced by Giidel 
et al. (2001a) to fit the flare observed during the first XMM- 
Newton observation of YYGem. For the modeling we fixed 
spectral components # 1 — 3 at their quiescent values (t\\ see 
Table Q). Generally, the statistics do not allow to constrain the 
abundances of the additional VMEKAL components. All abun- 
dances of components # 4 and # 5 were, therefore, held fixed at 
solar values. The 'high state' (£15) is an exception: Here the sig- 
nal at the high-energy end of the EPIC spectrum is better than 
for all other time-sliced spectra. Broad emission from the Fe K 
complex is clearly visible at ~6.7keV, and an acceptable solu- 
tion is obtained only if the Fe abundance is a free parameter, 
resulting in Fe/Fe = 0A7±%$. 

We point out that we performed an additional series of spec- 
tral fits where we set the abundances of these components to the 
quiescent values, and obtain much larger uncertainties for the 
model parameters. This hints at an abundance increase during 
active phases that seems to be not uncommon for stellar flares 



(Favata et al. 1999, Giidel et al. 1999) 



Our modeling allows to investigate how the temperature 
and emission measure of the flare components behave in dif- 
ferent stages of flares. In Fig. || we display both free tempera- 
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tures and emission measures from the 5-T model as a function 
of time together with the EPIC pn lightcurve and the evolution 
of the luminosity as derived from the spectral fits. Both spec- 
tral components show their highest temperature during the rise 
phase of the lightcurve. Subsequently, the temperatures decline 
coherently with the intensity. The emission measure increases 
during the rise phase and reaches its maximum with the begin- 
ning of the decay of the lightcurve. The observed luminosity 
L x ranges from (2 — 7.5) 10 29 erg/s in the EPIC band which 
is comparable to the values reported in the literature. The spec- 
tral parameters during extended time intervals with enhanced 
count rate (in, £15) suggest that these are flare-like events with 
a strong low-T component (VMEKAL #4) and a weaker high- 
T component (VMEKAL #5). But note that no typical flare 
outburst is seen in the lightcurve for the respective time inter- 
vals. 

5. Loop modeling 

The size of emitting coronal structures can be estimated by help 
of loop models. Most of these models are based on the assump- 
tion that the plasma is confined in a single loop. While this may 
be adequate during simple flares, the situation in the quiescent 
corona is likely to be much more complex. 

We used the temperature - emission measure approach 
which allows for heating during the flare decay, and determined 
the size of flaring loops from the spectral parameters of the 
EPIC spectra. We restricted the analysis to the first flare which 
is larger and longer than the second flare, and, therefore, pro- 
vides better sensitivity for this kind of study. 

The evolution of temperature and emission measure puts 
important constraints on the dynamics during flare decays. A 
one-dimensional hydrodynamic model describing the decay 



phase of solar flares was developed by Peres et al. (1987) and 



discussed in a series of papers (see e.g., Serio et al. 1991 



Jakimiec et al. 1992j |Sylwester et al. 1993 



Reale et al. 1997) 



The model was adapted to stellar flares by fceale et al. (1997' . 
The duration of the heating (a free parameter in the model) 
determines the slope £ in the lg T — lg V EM— diagram. The 
limiting cases are (a) abrupt switch off, and (b) a quasi steady- 
state evolution for slow decay of the heating which leads to 
evolution along the path given by the RTV scaling law (Rosner 
et al. 1978). The slope £ and the decay time of the lightcurve 
tlc are linked by the thermodynamic decay time scale, Ttd, i-e. 
the decay time if heating were absent. This relation depends 
on instrumental properties. With the appropriate calibration £ 
can be used to obtain an estimate for th e half-length L of th e 
flare loop, since according to Eq. (2) in Reale et al. (1997) L 
is connected to the maximum flare temperature, T max , and the 
thermodynamic decay time, r t d- 

The lg T — lg V EM -diagram derived from our modeling 
of the XMM-Newton EPIC spectrum during the larger flare is 
shown in Fig. [| As the model was developed for a single- 
loop and one-temperature flare, we ignored the first compo- 
nent (VMEKAL # 4) of the flare spectrum. This component 
has much smaller emission measure, i.e. the contribution to the 
flare emission is negligible compared to VMEKAL # 5. The 
curve in Fig. ^ starts with the rise phase (t 3 ) of the flare. We 



7.6 



7.5 



7.4 




25.7 



25.8 



25.9 
0.5 lg 



26.0 

'EM,) 



26.1 



26.2 



26.3 



Fig. 9. Temperature - emission measure - diagram for 
VMEKAL component # 5 for the larger of the two flares in 
the Sep 2000 XMM-Newton EPIC observation of YY Gem. The 
numbers next to the data points indicate the respective time in- 
tervals. 



measured the slope during the decay phase (t4...to) and found: 
C = 0.21. 

We applied the model described above. Modeling the de- 
caying part of the lightcurve by an exponential we found 
tlc = 16 ± 1 min. The relation between £ and rtd has 
been calibrated for EPIC-pn (Reale, priv. comm) and yields 
L = 2.010 9 cm for the length of the flare loop. As the low 
value we found for the slope £ hints at continuous heating, 
the evolution of the plasma parameters can in approximation 
be thought of as quasi-static. Then, using the number for the 
loop length given above, the RTV scaling laws imply a pres- 
sure of p ~ 6 10 3 dyne cm~ 2 for a hydrostatic loop, and a 
heating rate of En ~ 45 erg cm -3 s _1 . This yields a density 
in the loop of n ~ 6 10 11 cm" 3 . Assuming pressure equilib- 
rium the density at temperatures typical for the formation of 
O VII are estimated to be ^10 times as high. This contrasts 
with the low-density limit found from the analysis of the time- 
averaged O VII triplet, but is consistent with the high R ratio 
found during the flare. Therefore, the loop model underlines 
the indications for increased density during flares described in 
Sect. I 



6. Summary and Conclusions 

In September 2000 YY Gem was observed simultaneously by 
XMM-Newton and Chandra providing a possibility to cross- 
check the performance of the instruments on both satellites. 
We found that the line centers and the line fluxes measured by 
LETGS and RGS are in good agreement with only few excep- 
tions. One of these is the underestimation of the O VII reso- 
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Fig. 8. Time evolution of spectral parameters during the two flares on Sep 29/30, 2000. Shown are T4, T5, EM4, and EM^ 
from 5-T model fits to the joint EPIC pn + MOS spectra of YYGem together with the pn lightcurve (top panel) and the X-ray 
luminosity L x (bottom panel). The remaining three spectral components were held fixed at their quiescent values. Error bars are 
90 % confidence levels. 



nance line flux due to an absorption edge in the RGS which 
was not yet taken account of in the response matrix at the time 
of analysis (len Herder et al. 2002): For the resonance line of 
the O VII triplet the RGS flux is found to be ~ 70% of the 
LETGS flux, but the measurements of the two instruments are 
compatible with each other at the 3 a level. We note that similar 
differences have been observed between LETGS and RGS for 



Procyon (see |Raassen et al. 2002| ). A deviation of the flux for 
the Ne X Lya line between the two sides of the LETGS seems 
to point at a problem on the right side of this detector at the 
respective wavelength (12.13 A). 

The high-resolution spectrum of YY Gem is dominated by 
oxygen and neon lines and lines from intermediate ionization 
stages of iron, generally Fe XVIII and below. Only few lines 
from highly ionized iron (Fe XIX and higher) are seen at the 



long wavelength end of the LETGS spectrum. The characteris- 
tic FeXXI line at 128.73 A is detected but very weak, under- 
lining that the corona of YY Gem is not dominated by the high 
temperatures (> 4 MK) typically observed on RS CVn binaries 
(see e.g., ^lewe et al. 2001[ , |Audardet al. 200 1[ ). 

We have probed the optical depth in the corona of YY Gem 
using characteristic ratios of FeXVII line fluxes. Comparing 
these ratios with measurements for other stars known to exhibit 
very different activity levels shows that optical depth effects do 
not play a role in stellar coronae of all types. 

The correlation between X-ray luminosity and temperature 
in stellar coronae is well known (see e.g., Giidel et al. 1997). 



Now, the dispersive instruments of Chandra and XMM-Newton 
allow to confirm this relation using intensity ratios of individual 



emission lines. Ness et al. (2002b) show in their comparison of 
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Chandra observations of ten stars that YY Gem fits well into 
the relation between L x and temperature probed by the O VIII 
to O VII flux ratio representing a moderately active star. 

The only element in the spectrum of YY Gem with a suf- 
ficiently strong and unblended helium-like triplet is oxygen. 
The R = f/i value derived from the observation presented 
here is compatible with an earlier measurement by Giidel et al. 
(2001a). For the time-averaged spectrum the ratios of the oxy- 
gen triplet yield a temperature of T = 2 — 3 MK and an up- 
per limit to the electron density of n c < 2 10 10 cm -3 . During 
the time of observation YY Gem showed strong variability in- 
cluding two flares. The time-resolved analysis of the high- 
resolution spectrum reveals variations of R which are sugges- 
tive of an increase in density during times of enhanced count 
rate. Our measurement is one of the first direct indications 
for the increase of the coronal density during flares on stars 
other than the Sun derived from high-resolution X-ray spec- 
troscopy. A density measurement in the corona of AB Dor us- 
ing an XMM-Newton measurement of oxygen triplet lines, e.g., 
resulted in comparable dens ities during flare and during qui- 
escence ( Giidel et al. 2001b ). A detailed investigation of line 
ratios in the short time intervals representative for the different 
activity levels of YYGem is constricted by low S/N. The time- 
resolved high-resolution spectroscopy, furthermore, shows that 
during flares the absolute amount of (oxygen emission by) rela- 
tively low-temperature plasma increases. Lines from elements 
probing higher temperatures more typical for flares were not 
strong enough for this kind of analysis. 

The presence of higher temperatures in the corona of 
YYGem is, however, evident from our time-resolved spec- 
troscopy with EPIC. The temperatures for the global fit of 
the quiescent EPIC spectrum (2.4 MK, 7.4 MK, and 20.8 MK) 
comprise the temperature range measured from the analysis 
of oxygen lines, but go well beyond these values. Two ad- 
ditional VMEKAL components are necessary to describe the 
EPIC spectrum during flares. The corresponding temperatures 
and emission measures are highest during the rise phase. An 
attempt to reproduce the EPIC spectrum with a number of 
Gaussian profiles usin g the information fro m the RGS spec- 
trum was described by Haberl et al. (2002) . But the low reso- 
lution of EPIC makes this instrument insensitive to changes in 
the relative strength of nearby emission lines during different 
phases of activity. 

With plasma parameters derived from spectral modeling the 
size of the (flare) emitting str uctures can be inf erred. We used a 
hydrodynamic approach (see Reale et al. 1997 ) involving heat- 
ing during the decay phase to determine the length of flaring 
loops on YY Gem. The overall size of quiescent stellar coronae 
can be estimated making use of plasma density and tempera- 
ture derived from the analysis of high-resolution spectra (see 
e.g., Ness et al. 2001, Mewe et al. 2001, and Raassen et al. 



ing results from global fitting with EPIC. We find that loops 
are small (2 10° cm), on the order of a few percent of the stel- 
lar radius (i?* ~ 0.6 Rq) of each of the components in the 
YY Gem binary. This is consistent with the compactness of the 
corona established from eclipse mapping for the earlier XMM- 
Newton observation of YYGem (Giidel et al. 2001a). Small 



2002). However, loop lengths derived from individual lines re- 
fer only to a very limited temperature range and the relevance 
of the single loop model for the quiescent corona is question- 
able. In our case the spectral parameters derived from individ- 
ual lines are not useful for this kind of analysis, since the upper 
limits on the density translate to lower limits of the loop size 
according to the RTV-laws. Instead, we derived loop sizes us- 



loops point at magnetic activity on the individual stars rather 
than large scale structures associated with overlapping magne- 
tospheres. Furthermore, the RTV-laws point at high densities 
in the loop, confirming the evidences for low R— ratio during 
times of flares. 
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